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Abstract 

Understanding the dynamics of particles in turbulent flow is important in many environmental 
and industrial applications. In this paper, the statistics of particle pair orientation is numerically 
studied in homogeneous isotropic turbulent flow, with Taylor microscale Reynolds number of 
300. It is shown that the Kolmogorov scaling fails to predict the observed probability density 
functions (PDFs) of the pair rotation rate and the higher order moments accurately. Therefore, 
a multifractal formalism is derived in order to include the intermittent behavior that is neglected 
in the Kolmogorov picture. The PDFs of finding the pairs at a given angular velocity for 
small relative separations reveals extreme events with stretched tails and high kurtosis values. 
Additionally, The PDFs are found to be less intermittent and follow a complementary error 
function distribution for larger separations. 
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1. Introduction 


Understanding how particles are advected by fluids is of major interest in many applications, 
including the environmental and the geophysical flows [l| . One outstanding example is the erup¬ 
tion of volcanoes, where particles with different sizes and inertia are released into the atmosphere 
and then transported with turbulent currents [^, |^. In addition to the environmental impor¬ 
tance of understanding the influence of turbulence on particle dispersion, this problem is_ also of 
practical interest in the industrial flows, where the advection of particles is involved 


The theoretical studies o 
models. Many experiments 


the relative separation between two particles are based on stochastic 
3,1^ and numerical simulations have been performed over the 

ast few decades in order to accurately evaluate these theories. Richardson, in his pioneering work 
Hi, has examined the relative motion of particle pairs in turbulent flows. He has estimated the 
scale dependency of the eddy-diffusivity coefficient through the observation of dispersed plumes. 
This scale dependency is claimed to be the origin of the accelerated nature of turbulent dispersion 
[l3 | . In the inertial range of motion (for r/ <C r <C T and <C t Tl , where rj is the Kolmogorov 
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dissipative scale, L is the energy injection length scale of the flow, is the local-eddy-turn-over- 
time and Tl \s the large time scale), Richardson has suggested that the process of relative 
dispersion is governed by a diffusion-like equation. Solving this equation gives the probability of 
hnding a pair of particles at a specific separation, at any time [1^. Nevertheless, recent studies 
|14l . Il5( on the dynamics of tracer pairs that are released from many point sources have reported 
severe deviations from the Richardson theory. Additionally, this theory does not account for the 
rotational behavior of particle pairs. 

The orientation dynamics of a single rigid ellipsoid particle under Stokes flow has been de¬ 
scribed by Jeffery 16|], where inertia and the thermal fluctuations are neglected. The same 
equation can be used to describe the motion of any axisymmetric particle, provided that its 
aspect ratio is known 17|. In the presence of weak inertia, Einarsson et al. 1^, ll9| studied 


the rotation of small and neutrally buoyant axisymmetric particles in a viscous shear flow, by 
perturbatively solving the coupled particle-flow equations. Shin and Koch [20| presented the re¬ 
sults of direct numerical simulations (DNS) of the translational and rotational motions of fibers 
in a fully developed isotropic turbulent flow, for a range of Reynolds numbers. They concluded 
that the fibers whose lengths are much smaller than the Kolmogorov length scale r/ translate like 
fluid particles and rotate like material lines. With increasing fiber length, the translational and 
rotational motions of the fibers slow down as they become insensitive to the smaller-scale eddies. 

Using computer simulations, Pumir and Wilkinson 2]| studied the temporal evolution of 
the orientation vector of microscopic rod-like particles. They showed that rod-like particles 
are more strongly aligned with vorticity than with the principal strain axis. An interesting 
component of turbulence, which is named the pirouette effect, has been reported by Xu et al. 


22l |. It has been shown that the axis of rotation of tetrahedra tracers aligns with the initially 


strongest stretching direction, a phenomena which can be justified by the conservation of the 
angular momentum. Using video particle tracking technique, Parsa et al. 23| made accurate 
simultaneous measurements of the motion and the orientation of the rods in the presence of 
a bi-dimensional chaotic velocity field. The first three-dimensional experimental measurements 
of the orientation dynamics of rod-like particles was also reported by Parsa et al. 2^. In 


this work, a good agreement with the previous numerical simulations was reported. Moreover, 
Parsa and colleagues 2^ studied the rotation rate of the rods with lengths in the inertial range 
of turbulence. They presented the experimental measurements of the rotational statistics of 
neutrally buoyant rods, and derived a scaling law for the mean-squared rotation rate. The latter 
showed a good agreement with the Kolmogorov classical scaling. 

It is noteworthy that, to the best of our knowledge, the orientation of the tracer pairs has not 
been studied so far. Therefore, this question will be addressed here, by computing the probability 
of finding a pair of particles at a specihc rotation rate, given the relative separation between them. 
By considering a particle pair with a specihc separation, one can study the orientation statistics 
of an imaginary tracer rod that is delimited by these two particles. Furthermore, for the hrst 
time, the pair rotation rate PDFs at later times will be presented in this work. 
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This paper is organized as it follows. The next section describes the numerical data, and the 
approach that is used to derive the rotation rate of the tracer pairs. Afterward, the intermit- 
tency behavior that is observed at the small scales of motion is discussed in detail. In Sec. [3l 
the multifractal (MF) formalism to intercept the probability density function (PDF) distribu¬ 
tions is presented. Thereafter, the higher-order moments of the rotation rate are evaluated and 
compared to the multifractal prediction together with the Kolmogorov (K41) picture. Finally, 
the concluding remarks, and the outlooks are stated in Sec. HI 


2. Methods 


2.1. Particle pair orientation 

The rotational dynamics of particle pairs in turbulent flow can be addressed by solving the 
dimensionless incompressible Navier-Stokes equations for a Newtonian fluid, i.e. 


^U. 1 

^-hu.Vu = -VpP —'V'^u + i, 
ot Re 

V.u = 0, 


( 2 . 1 ) 

( 2 . 2 ) 


using Direct Numerical Simulations (DNS). In these equations, u denotes the tri-dimensional 
velocity field, p is the pressure and Re is the flow Reynolds number, Re = Lv!/v. Re measures 
the ratio between the nonlinear inertial forces, and the linear viscous forces, where u' is the 
root-mean-squared velocity and v is the fluid kinematic viscosity which is defined as the ratio 
between the dynamic molecular viscosity /U and, the fluid density p. In order to avoid energy 
dissipation, the flow was forced by keeping the total energy constant in the first wavenumber 
shells, by applying a large-scale forcing term f. This force injects ener gy at a mean rate of 
e = (f.u), where e is the mean turbulent kinetic energy dissipation rate |26l |. The integration 
of the equations of motion is performed on 1,024^ regular cubic lattice with Taylor microscale 
Reynolds number of 300, and periodic boundary conditions. The simulation parameters are 
shown in Tab. [TJ With the present choice of parameters, the dissipative range of length scales 
is well resolved because the grid size Ax is in the range of Kolmogorov length scale rj (as it 
is reported in more detail by Biferale et al. (3). A fully-de-aliased parallel pseudospectral 
code, with a second-order Adams-Bashforth temporal scheme, for 3D homogeneous isotropic 
turbulence, assuming constant fluid viscosity and density, is used in order to solve Eqs. m 
and ()2.2p . 

The fluid is seeded with bunches of tracers emitted within a small region which has a size 
comparable to the Kolmogorov dissipative scale p. The emission is carried out in puffs of 2,000 
particles that are followed during a maximum time of 160 Tr^. The tracer particles take on the 
fluid velocity immediately, and adapt the rapid fluid velocity fluctuations [27|. Therefore, the 
velocity of each tracer is related to the instantaneous fluid velocity by the following equation 


dxr 


Vp ^ 


dt 


= u(xp(t),t). 


(2.3) 
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Table 1: Parameters of the numerical simulation. 


Kolmogorov dissipative active length scale r] 

0.005 

Spacing between two collocation points in the regular cubic lattice Ax 

0.006 

Mean turbulent kinetic energy dissipation rate per unit mass e 

0.81 

Fluid kinematic viscosity v 

0.00088 

Dissipative time scale 

0.033 

Integral time scale Tl 

67 

Energy injection length L 

27r 

Number of collocation points N 

1,024 

Turbulent velocity fluctuation u' 

1.7 


The particle trajectories are computed by integrating the Eq. (|2.3I) . The positions and the 
velocities of each particle are stored at a sampling rate of For the 2,000 particles generated 
in each puff, all the possible pairs (approximately two millions) are considered 281. 


2.2. Pair angular velocity 

The pair angular velocity can be measured on the basis of the instantaneous positions and 
the velocities of the two particles. Consider a given pair of particles defined at every time by 
its separation vector r pointing from the first tracer particle toward the second one. One of the 
particles is taken as a reference point for computing the pair orientation. The relative velocity 
of the other particle Au is decomposed into two components; i.e. Au|| parallel to the separation 
vector, and Auj^ perpendicular to it. The first particle and the transverse component of the 
relative velocity of the second one, defines a plane of rotation. The axis of rotation e is then 
normal to this plane, and dehnes the direction of the angular velocity pseudovector cj. By taking 
9 as the angle between the separation vector r and Au, then the angular velocity vector can be 
written as 


u> = 


or 


and in tensorial notation. 


Auj_ Au sin0 

—i—i— e =-—j- e. 

(2.4) 

r r 

r X Au 

(2.5) 

U) = - 

^2 

2 ’ 

(2.6) 


where Cijk is the Levi-Civita symbol. Considering the isotropy of turbulence, the ensemble 
average over the angular velocity components vanishes, that is to write (ujk) = 0. 


3. Results and discussions 
3.1. Multifractal model 

The angular velocity of the rotating particle pairs in Eq. (|2.6p can be computed at any 
time by knowing the spatial positions, and the velocities of the two tracer particles separated 
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by a vector r. In the case of inertialess particles, the angular velocity can be viewed as the 
ratio between the Eulerian transverse velocity increment 6rU± = u±{r + 6r) — u±{r), and the 
instantaneous distance r. The mean-squared angular velocity (w^) for a fixed separation r, is 
thus the ratio between the second order Eulerian transverse velocity structure function, and 
Eor the purpose of computing the higher-order moment of angular velocity that is parametrized 
by r, a series of computations over the whole number of bunches at different times, and for dif¬ 
ferent constant separations are performed. In the inertial range of motion, the size of the eddies 
that move the pairs apart varies with the separation distance. The eddies with the scales in the 
order of r are the most effective in the dispersion process 29|]. Intuitively, in this intermediate 
range, the further the separation between the two particles is, the slower the pair will rotate. 
This happens because in the Kolmogorov phenomenology, the angular velocity modulus scales as 
j,-2/3_ According to the Kolmogorov second similarity hypothesis 30], all statistically averaged 
quantities at scale r in the intermediate range are uniquely determined by e and r. The pth 
order moment angular velocity should then scale in the inertial range as 


{ujP) ~ r“2p/3_ 


(3.1) 


However, in the dissipative range of turbulence (i.e. for r rj), the higher order moments 
are independent of the scale r. A suitable expression is therefore needed in order to fulfill the 
requirements of both ranges. In Kolmogorov phenomenology, the higher order moments are 
expressed with respect to separation r via the Batchelor parametrization which satisfies 
both the dissipative range and the inertial range. Then (ujP) is expected to be scaled as 32] 


(ujP) ~ 


(x) 


-2p/3 


1+1^ 


-2\ P/3 ’ 


(3.2) 


Eor separations in the inertial dissipative range (i.e. r > r] and r L), these moments 
undergo a quadratic decay. 

Kolmogorov theory fails to accurately predict the intermittency behavior, because it neglects 
the presence of fluctuations in the energy transfer from the large scales toward the small ones 
33|. In order to include the intermittency effect in the energy cascade, a multifractal formalism 
is used together with a generalization of the Batchelor parametrization, to describe the higher 
order moments of the angular velocity for both the inertial dissipative range and the inertial range 
34l |. The multifractal approach to turbulence is based on the assumption that the statistical 


properties of turbulent flows do exhibit scaling properties even i 
Under these assumptions, the transverse velocity increment reads 


there is intermittency 35l |. 


32] 


6rU± = Vo 


[iY 




_2X {l-h)/2' 


(3.3) 
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where fo denotes the large scale transverse velocity fluctuation and h is the singularity in the 
MF phenomenology. In this picture, the viscous scale is not a fixed homogeneous quantity, but 
it fluctuates wildly due to intermittency 3^ . The Kolmogorov length scale r/ is a function of the 
singularity h, and it is related to the Reynolds number by 77 /L oc 

Following Eqs. (j2.6p and (I3.3p . the angular velocity component at a distance r is estimated 
in the MF model as 

,h-l 

(3.4) 


OJ = 


6rU± 


L ' 


ay 




_2X (l-/*)/2’ 


A generalization of the probability to observe an Eulerian velocity fluctuation at scale r with 
a local Holder exponent h, is given by 3^ 


Ph{r) (X 


ii) 


3-D{h) 


1 +' '■ 




_2X (D(h)-3)/2’ 


(3.5) 


where D{h) is the spectrum of the fractal dimension. 


3.2. Probability density function distribution 

In this study, there are many particle pairs which are embedded in the turbulent flow. The 
most comprehensive and accurate description of the orientation state can be achieved by evalu¬ 
ating the PDEs of the orientation of the pairs. The non-Gaussianity is one of the fundamental 
features of 3D turbulence with a non-zero skewness, and a very large kurtosis of the PDE distri¬ 
butions (37|. A superposition of stretched exponentials corresponding to the various singularity 


exponents are usually observed at small scales. In the following, we shall be interested in the 
PDE of hnding a pair of tracers at a given angular velocity component oj when it reaches a given 
separation r. 

The ME phenomenology is capable of predicting the PDEs of the velocity increment and 
the accelerations in both the Eulerian and Lagrangian frames 


38 


39|. The ME model can also 


successfully intercept the numerically observed angular velocity PDEs. In this section, an ME 
model is employed for the PDEs of the angular velocity components. The PDEs are allowed 
to take values in all the axis of reals. Since the flow is isotropic, i.e. the pair properties are 
distributed independently of orientation, the distributions are symmetric. Eurthermore, the 
PDEs that correspond to each of the components all collapse on a single curve. 

By assuming that the large scale velocity field vq is Gaussian with d components, its magni¬ 
tude has the following PDE distribution 40l | 


..d-l 


P{vo)dvo = Vq " exp I -y 


(3.6) 


where the normalization constant has been omitted here. The transverse velocity modulus corre- 
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Figure 1: (Color online) A comparison between the numerical results (in gray) and the MF prediction of the 
normalized PDFs of the angular velocity component, that are calculated for the separation lags of r/r] = 1, 40, 
200 and 1,000. The dashed line represents the unit variance complementary error function given by (13.811 . It can 
seen that the complementary error function overlaps perfectly on the numerical results for r = 1, OOO 77 . Inset: The 
variations of the kurtosis of the angular velocity component PDF. The horizontal dashed line corresponds to 18/5 
as the kurtosis of the complementary error function. 


spends to d = 2, and the PDF follows a Rayleigh distribution. By setting the variance to unity, 
one can write the large scale transverse velocity PDF modulus as 

P{vo) = uoexp ^ (3-7) 

It is possible to derive the PDF of the large scale transverse velocity component fluctuation 
from the PDF of the transverse velocity modulus given by (13.7p . Because the three components 
have the same probability distribution due of isotropy, the PDF of each of them is given by a 
complementary error function (see the appendix for more details). For a variance of unity the 
following large scale transverse velocity component PDF is derived 

The probability of finding a pair of particles at a specific angular velocity, and for a given 
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separation r is obtained from 


P{oj\r) 


= f P{oj\h)Ph{r) dh. 

h£l 


(3.9) 


By plugging Eqs. (|3.8p and (I3.5p into (|3.9p . and using the probability conservation law for 
the equation (I3.4p . one can write 


P{Lj\r) oc 



D(h) + h 


erfc 


h£l 



' (r\ 

1 V3 1 

[lj 



(3.10) 


The fractal dimension has the following form in a log-Poisson process 

3(/i - ho) 


41 


m 


D{h) = 


ln/3 


■ / 3(ho - h) 


- 1 


+ 3 — do, 


(3.11) 


where /3 = 0.6, ho = 1/9 and do = (1 ~ 3/io)/(l — /3) = 5/3. The lowest singularity hmin = 1/9 is 
obtained by considering lim/j^/j^.^ D{h) = 0. The fractal dimension attains its maximum which 
is D — 3, for hmax — hp^ak — 0.38 36]. 

In Kolmogorov theory, the singularity h takes only the value of 1/3, and the fractal dimension 
is restricted to D = 3. It is possible to derive the PDF distribution in this particular case from 
Eq. ()3.10p . After computation, it is found that the angular velocity component PDF is given by 


P{uj\r) 


^^erfcf— 

^rms 2V3 VV 3 ^rms / 


(3.12) 


where Wrms is the root-mean-squared angular velocity, expressed in Eq. (13.2p for p = 2. The com¬ 
plementary error function distribution given by (|3.8h . is therefore recovered for the normalized 
PDF. 

3.2.1. Comparison with numerical results 

In Fig. [H a comparison between our numerical results and the MF prediction given by (|3.10l) 
is presented. The agreement between these two shows the robustness of the MF approach in 
reproducing the DNS results accurately. The PDF of the angular velocity component for the 
smallest separations i.e. of the order of the Kolmogorov length scale p, is highly intermittent and 
the PDF has stretched tails that extend beyond thirty times the standard deviation of angular 
velocity (for the sake of convenience, in Fig. [T]tails are not shown in full lengths). This indicates 
that when the separation length r is of the order of the smallest scale of the flow, the probability 
of high rotation rate events are much higher than those predicted by a complementary error 
function distribution. An increasingly large number of particles is needed in order to ensure a 
sufficient statistical convergence of the data, and therefore to describe the long tails accurately. 
The inset of Fig. [T] illustrates this by showing the variations of the PDF kurtosis as a function of 
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Figure 2: (Color online) Comparison between the Kolmogorov scaling and the multifractal model. The scaled 
higher order moments of the angular velocity component is depicted versus the scaled separations in the inertial 
range. For p = 2,4, the MF model and the K41 scaling given by are almost indistinguishable. The 

multifractal prediction shows a better agreement with the numerical results for the higher order moments p = 
6 , 8 , 10 , 12 . 


r. The kurtosis is a measure of the distribution pickedness, 
it is defined as the ratio between the fourth order moment, 
variance. That is to write, 


And for symmetrical distributions, 
and the square of the distribution 

(3.13) 


For small separations, the kurtosis is very large and it reaches a value of 10 for a separation 
of the order of the Kolmogorov scale. The larger the separation is, the smaller the kurtosis will 
be, and it finally converges toward 18/5. This is the kurtosis of the large scale angular velocity 
component PDF given by (13.8p (see the appendix for more details about the computation of this 
kurtosis). 


3.3. Higher order moments 

The pth order moment of the angular velocity can be computed from the PDF distributions 
as 

(w^) = J LJ^P{Lo\r)duj. (3-14) 

ojeM 

By plugging Eqs. (j3.10l) into (|3.14l) . and integrating uj over the whole axis of reals, the integral 


9 














expression of the higher order moment of the angular velocity component at separation r is 



hei 



|-t(D(/i)+p(l-/i)) 

dh. 


(3.15) 


When the support I is limited to h = 1/3 with a fractal dimension of II = 3, the Kolmogorov 
scaling given by (13.2h is recovered. Eq. (j3.14j) can also be written as follows, by integrating with 
respect to the singularity h 

{u}^) ~ f LoPP{u!\h)dh. (3.16) 

h€l 

In Fig. [2] the scaled higher order moments of the angular velocity component are plotted 
against the scaled separations r that extend from the inertial range deep into the dissipative 
inertial range. For the first two even order moments, the Kolmogorov picture and the MF 
prediction are almost indistinguishable. In fact, the intermittency behavior can only be observed 
for the higher order moments, for which the simulation results present a steeper scaling. 

Nevertheless, the Kolmogorov prediction given by (|3.2I) cannot accurately reproduce the 
numerical results for the higher moments beyond the sixth order moment. Additionally, the 
normalized PDFs that are computed using this theory, for different separations all overlap to 
the same graph. In other words, K41 can only be applied for the large separations, when the 
complementary error function distribution is adopted. In contrast, the multifractal model is 
found to be in a satisfactory agreement with the simulation results for all the separation ranges. 

Assuming that in the proposed MF model, the pih. order moment of the angular velocity in 
the inertial range, can be written as 


(cuO 


r.~fp + U) 


(3.17) 


that is to write the scaling exponent as the sum of the exponent which is predicted by the 
Kolmogorov theory in Eq. and a correction which takes into account the intermittency 

effects. For the higher order moments, the correction has a negative contribution in the exponent, 
and the MF model leads to a steeper scaling, giving a better agreement with the numerical results. 
This happens because MF includes the intermittency character of the energy dissipation which 
is neglected in the classical theory. 

In Fig. [3] the PDF distribution of the pair angular velocity components at various times 
is illustrated. At the initial time, the particles are located at a source with dimensions of the 
order of the Kolmogorov scale rj. Scatamacchia et al. 15|, 1431] reported an abrupt transition in 


the particle dispersion at a time of about lOr^ after the bunch emission. At that time, most of 
the pairs reach a relative distance of the order of lOr/. At times larger than 20rr;, the relative 
separation of the majority of pairs is within the inertial range of turbulence. In this range of 
motion, the angular velocity is scale dependent, and it is inversely proportional to the distance 
separating the two particles. This means that at later times, the probability of finding the pair 
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Figure 3: Lin-log plot of the positive part of the symmetrical PDF of finding the pair of tracers at a given scaled 
angular velocity, in the range extending from 20 (solid line) to flO r^. In this figure, to is normalized by 
the inverse of the eddy-turn-over time at the Kolmogorov scale r^. At later times, the bunches experience more 
expansion, and the probability of finding the pairs at a lower rotation rate becomes higher. 

at a smaller angular velocity is much higher, because the bunch of particles has experienced a 
larger dispersion. The maximum of the PDF distribution increases with time and the region of 
small rotation rates becomes more populated. 

4. Conclusion 

In this work, the orientation statistics of two particles as they are carried in a flow are studied. 
The flow is seeded by tracer particles that are released from one point-like source in bunches of 
thousands. Knowing the instantaneous positions and the velocities of the emitted particles, the 
higher order moments of the angular velocity component of the rotating pairs are measured. The 
theory of Kolmogorov based on purely dimensional arguments is found unable to accurately fit 
the numerical results. 

The multifractal phenomenology is therefore employed in order to include the intermittency 
corrections. It is shown that it could remarkably reproduce the observed behaviors for the 
angular velocity higher order moments and for the orientation PDFs. The computation of the 
PDFs of finding the particle pairs at a specific angular velocity component when they reach 
a given separation reveals that one can quite easily attain events in excess of thirty times the 
standard deviation with high kurtosis values. The dramatic change of the angular velocity PDF 
shape provides with immediate impression of how strongly intermittent rotation events are at 
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small scales of motion. By increasing the separations, the PDFs become less intermittent and 
converge toward the complementary error function distribution. 

Here, the main attention is restricted to the neutrally buoyant particles which are passively 
transported by the turbulent flow. The role of particle inertia on pair orientation should be 
explored in the future. It is important to note that the measurements in the database are not 
completely unbiased, because only the particles that are emitted from a single point source are 
considered. In order to reduce the spatial correlation, it is worth considering the pairs which 
are formed from various point-like sources. For reliable statistics, a large number of particles 
will be required. It is worth mentioning that an ideal framework for the study of the orientation 
statistics is a 2D turbulence, because of the absence of intermittency. Similar measurements and 
quantitative investigations can be carried out as well in this particular case, where the angular 
velocity vector has a unique component along the axis normal to the computational domain. 
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Appendix 


Here, the PDF of the components are derived for a given vectorial quantity, by knowing the 
PDF of the modulus when the isotropy is guaranteed. In the following, x, y and z represent the 


vector components. Let N{r) be the distribution function of the variable r where = x'^+y'^+z'^. 
N{r) is normalized to unity, that is to write 


Using the spherical coordinates, Eq. can be written as. 



(5.1) 



(5.2) 


The PDF for the modulus Q{r) = 47rr^A(r) has a total weight of unity. Because of isotropy, 
the PDFs are all the same along all the three axis. Then only P{z) will be calculated here. 



(5.3) 
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P{z) can be written as a function of Q{r) 


P{z) 


1 

47r 



+ 2/2 + ^ 2 ) 172 ^ 

_|_ y2 _|_ ^2 


dxdy, 


(5.4) 


If the two-dimensional vector whose Cartesian coordinates are x and y are described in the 
cylindrical coordinates (u, (j), z), where u varies over [0, oo[ and (j) over [0, 27r], then + 2/^) 

and 


1 r 


oo Q( [u^ + 2;^ 


udn, 


(5.5) 


v? + 2;2 

By considering the change of variables = v? + z"^, then tdt = udu, the following PDF for 
the positive values of z will be obtained, 


1 r 


Qit) 


dt. 


(5.6) 


The resulting equation (15.61) can be illustrated in one well-known example of turbulence. It 
is known that the large scale velocity has the following distribution, 


Q{t) = -^t^e * . 


TT 


(5.7) 


Using Eq. (15.61) one finds that the component PDF has a Gaussian distribution, 


P{z) = 

/tt 


(5.8) 


which has a total weight of unity. 

As it is discussed in 13.21 , the large scale transverse velocity fluctuation has a Rayleigh distri¬ 
bution given by Eq. (13.7p . then the PDF of each component has a complementary error function 
distribution after integrating Eq. (15.61) . 


p(-) = 


^ertc(' 


Z\ 

\V2 


(5.9) 


which has also a total weight of unity. 

Moreover, using Eq. (|3.13p the kurtosis of the above distribution is calculated as 


K = 




^2 erfc ( 

VV2 


Y' 


(5.10) 
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